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1. INTRODUCTION 


The Agricultural and Resources Inventory Surveys Through Aerospace Remote 
Sensing (AgRISTARS) is a 6-year program of research, development, evaluation, 
and application of aerospace remote sensing for agricultural resources begin- 
ning in fiscal year (FY) 1980. The AgRISTARS program is a cooperative effort 
of the National Aeronautics and Space Administration (NASA), the U.S. Agency 
for International Development (AID), and the U.S. Departments of Agriculture, 
Commerce, and the Interior (USDA, USDC, and USDI). 

The goal of the program is to determine the usefulness, cost, and extent to 
which aerospace remote sensing data can be integrated into existing or future 
USDA systems to improve the objectivity, reliability, timeliness, and adequacy 
of information required to carry out USDA missions. The overall approach is 
comprised of a balanced program of remote sensing research, development, and 
testing which addresses domestic resource management as well as commodity 
production information needs. 

The technical program is structured into eight major projects as follows: 

1. Early Warning/Crop Condition Assessment (EW/CCA) 

2. Foreign Commodity Production Forecasting (FCPF) 

3. Yield Model Development (YMD) 

4. Supporting Research (SR) 

5. Soil Moisture (SM) 

6. Domestic Crops and Land Cover (DCLC) 

7. Renewable Resources Inventory (RRI) 

8. Conservation and Pollution (C/P) 

The majority of these projects will make direct use of information on crop 
phenology. Specific areas of these projects to which phonological information 
is pertinent include classification, acreage and yield estimation, and detec- 
tion of episodal events. 
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This report presents results of developmental studies directed tov/ard the: 
identification of crop calendar techniques suitable for use in AgRISTARS 
projects. Although most of the approaches are generic in nature, the focus 
has been on spring wheat, with the expectation that the initial user in the 
AgRISTARS program will be the FCPF project. 



2. BACKGROUND 


Prediction of crop growth stages (phenology) may be achieved through three 
independent methods: calculating historical averages for an area (normal crop 

calendars); agrometeorological modeling from knowledge of crop-weather inter- 
actions [agricultural-meteorological (agromet) crop calendars]; and interpre- 
tation of spectral signatures (spectral crop calendars). Models of these 
types have been and are currently being developed. The quantification of 
crop-weather-spectral relationships (spectromet crop calendars) in a single 
model is still in the early stages of research. 

The nature of the crop stage estimation problem, however, allows for the 
development and testing of crop calendar systems which have components of dif- 
ferent types. Typically, an agromet crop stage model requires knowledge of 
planting dates. Thus, a spectral approach which estimates planting dates 
could be used to initialize an agromet stage model. This kind of system may 
be considered a first generation spectromet model. 

The following paragraqihs provide an overview of the general characteristics of 
the major crop calendar components and systems which were examined in this 
study. 

The normal crop calendars provide a baseline against which other methodologies 
can be evaluated. A simple measure of skill lies in the amount to which a 
candidate technique reduces the variance in stage estimates as compared to the 
reduction achieved by normal crop calendars. 

Adjusted normal models utilize planting date estimates to shift the normal 
crop calendars in time. Agromet planting date models providing segment-level 
estimates were evaluated both independently and as part of such hybrid systems. 
Due to constraints in the available data, spectral planting date estimation 
systems providing field-level estimates were developed and evaluated only in 
the context of the adjusted normal model. 
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Normal, agromet, and spectral planting date models were evaluated in conjunc- 
tion with two agromet crop development models due to lack of data sets 
suitable for independent evaluations. 

Two spectral techniques which made direct estimates of crop growth stages on 
each Landsat acquisition date were developed and evaluated. The models were 
evaluated on a subset of 1978 crop-year data for which there were simultaneous 
ground truth and spectral observations. Each acquisition was considered 
independently for both techniques. 

The spectral stage estimates for each acquisition were also used with a normal 
crop development scale to estimate planting dates. The average spectral 
planting date with respect to field was calculated for use in the previously 
mentioned spectromet and adjusted-normal models. 

In contrast, but also falling into the adjusted- normal classification, a spec- 
tral technique which fit the multitemporal signature sequence was evaluated. 
For this technique, at least three spectral acquisitions for application were 
required. The technique provided an estimate of the date of maximum spectral 
response. The date of peak response was combined with a normal crop develop- 
ment scale to estimate the planting stage date as well as other growth stage 
dates. This approach was evaluated against ground truth observations of 
growth stages for those segments in which it could be applied. 

The following sections provide a summary of the data used, specific models 
used, and evaluation results. In most cases, technical details of the models 
and developmental results have been placed in appendices of this document. 
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3. DATA 


Two distinct data sets were available for use in evaluating crop stage models. 
The first set was collected at intensive test sites (ITS) of the Large Area 
Crop Inventory Experiment (LACIE) during crop years 1975 to 1977 (ref. 1). 

The second set was collected at the LACIE Transition Year (TY) blind sites 
during the 1978 crop year (ref. 2). 

Limitations inherent in each of these data sets placed constraints on the 
nature of possible test runs. The data needed included field-level observa- 
tions of planting date, periodic observations of phenology during the season, 
and multitemporal Landsat observations of the same fields along with repre- 
sentative daily meteorological data. The ITS data set was found to be 
suitable only for testing agroment planting date models. The number of 
spectral acquisitions was too small for evaluation of spectral models. The 
phenol ogi cal observations were also few in number and were taken on a scale of 
limited usefulness. The ITS data set is shown in table 3-1. 

The TY data set was complete except for inclusion of the planting dates. This 
data set was used to evaluate the spectral stage estimation models and the 
hybrid model systems. Data were collected for spring wheat at blind sites in 
Montana, North Dakota, South Dakota, and Minnesota. Ground truth observations 
were collected for special fields every 9 or 18 days. The development of the 
crops was observed using the Feekes scale. Table 3-2 gives a summary of the 
number of fields observed, the number of usable Landsat acquisitions, the 
number of ground truth observations, and the number of fields available for 
evaluation with the multitemporal profile (Badhwar) technique. The locations 
of the segments are shown in appendix A. 
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TABLE 3-1.- COMPONENTS OF EVALUATION DATA SET 
FOR PLANTING DATE MODELS 


State 

Segment 

Year 

Number of 
fields 

North Dakota 

1965 


8 

Minnesota 

1987 


45 

South Dakota 

1687 

1976 

9 

North Dakota 

1965 

1976 

. 8 

North Dakota 

1966 

1976 

31 

North Dakota 

1967 

1976 

5 

Montana 

1970 

1976 

2 

Mi nnesota 

1987 

1976 

12 

South Dakota 

1687 

1977 

10 

Minnesota 

1987 

1977 

11 
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TABLE 3-2.- BLIND SITE TRAIN ING/EVALUATIOM SET FOR SPRING WHEAT 


State 

Segment 

Number of 
fields^ 

Number of 
Landsat 
acquisitions 

Number of 
periodic 
observations 

Number of 
usable 

observations 

N. Dak. 


, 5 

2 

5 

25 

N. Dak. 


^11 

3 

5 

55 

N. Dak. 


8 (3) 

6 

4 

29 

N . Oak . 

1457 

5 (5) 

5 

6 

31 

N. Dak. 

1461 

13 (12) 

7 

8 

70 

N. Dak. 

1472 

2 

3 

7 

13 

N. Dak. 

1473 

4 

3 

4 

16 

S. Oak. 

1485 

13 

3 

3 

41 

Minn. 

1518 

12 (12) 

7 

8 

79 

Mont . 

1537 

2 (2) 

7 

6 

11 

Mont . 

. 1542 

6 (6) 

6 

5 

29 

Mont. 

1544 

4 (3) 

5 

5 

18 , 

Mont. 

1553 

3 (3) 

7 

7 

18 

Minn. 

1566 

15 (8) 

4 

7 

97 

N. Dak. 

1584 

8 

3 

6 

40 

N. Dak. 

1612 

7 (3) 

4 

6 

39 

N. Oak. 

1619 

8 (7) 

4 

5 

40 

N. Dak. 

1636 

10 (10) 

9 

7 

54 

N. Dak. 

1658 

7 

2 

6 

40 

N. Oak. 

1664 


- 3 

6 ’ ' 

■ - - ^ 35 - - 

S . Oak . 

1668 

7 (3) 

4 

6 

38 

S. Dak. 

1755 

3 (2) 

7 

5 

14 

S. Oak. 

1811 

12 (10) 

5 

6 

58 

Minn. 

1825 

13 (13) 

5 

4 

52 

N. Dak. 

1918 

4 

3 

6 

21 

N. Oak. 

1920 

8 (7) 

6 

6 

47 

N. Oak. 

1924 

4 (4) 

8 

6 

22 

Mont. 

1942 

’6 

2 

5 

25 

Total 

28 

206 (113) 



1057 


^Numbers In parentheses indicate number of fields usable for the Badhwar 
multi temporal profile fit. 
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4. MODELS 


Three estimators of spring wheat planting dates were evaluated using as ground 
truth the planting dates of the fields in the ITS data set. One was the 
normal planting date as estimated from historical Crop Reporting District 
(CRD) dates. For Minnesota and South Dakota, the estimates are based on ten 
years of data. Only five years of data were available for North Dakota and 
Montana. 

Two agromet planting date models, the Feyerherm and the Stuff-Phinney, were 
evaluated. Details of these models may be found in appendices B and C, 
respectively. Using weather data to estimate the rate of planting, segment- 
level planting dates were predicted by these models. 

The model systems tested are described in the following paragraphs 
(sections 4.1 through 4.7). The available data set did not contain observa- 
tions of planting dates. Therefore, a planting date estimation algorithm was 
tested in conjunction with a crop stage estimation algorithm. The mnemonic 
labels for each system indicate the starter model and the crop stage model that 
were used. 

4.1 THE NORMAL-NORMAL MODEL 

The basic concept of this model is that historical averages can be used to 
predict current year crop development. The mechanism by which this is done 
consists of applying the standardized crop development scale, described in 
appendix D, to the historical averages for the CRD's in which the various seg- 
ments are located. The normal planting and harvest dates applicable to each 
segment are given in appendix E. 

4.2 FEYERHERM-AD JUSTED NORMAL MODEL 

The basic concept of this model is that the Feyerherm agrometeorological 
planting date model can be used to shift the normal model's standardized crop 
development scale (appendix D) backward or forward in time and thus account 
for a major portion of the year-to-year variability in crop development rate. 
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4.3 FEYERHERM-ADJUSTED ROBERTSON MODEL 


The basic concept of this model is that the combination of two agrometeorolog- 
ical models will provide better estimates of growth stages than would either 
one used in conjunction with the normal model. The Robertson model is a 
triquadratic multiplicative model that uses photoperiod and maximum and 
minimum daily temperatures to estimate growth stages. It is initiated on the 
planting date estimated by the Feyerherm model. A description of the 
Robertson model is presented in appendix F. 

4.4 FEYERHERM-ADJUSTED DORAISWAMY MODEL 

The basic concept of this model is similar to the preceding one, with the 
additional assumption that performance of the Robertson model can be improved 
by making the coefficients more realistic than those originally derived by 
statistical fitting. Details on the Ooraiswany agromet model may be found in 
appendix G. 

4.5 BADHWAR- ADJUSTED NORMAL MODEL 

The basic concept of this model is that the model function, when fit to field 
spectral data, provides peak Kauth greenness dates which will improve the per- 
formance of the normal model. Further details on the development and adapta- 
tion of this model are provided in appendix H. 

4.6 COLOR-ADJUSTED NORMAL MODEL 

The basic concept of this model is that an empirical logic based on color- 
coded Landsat data can be combined with the normal model to provide improved 
estimates of growth stages. The independent performance of the color logic in 
predicting coincident ground truth biostages was also measured. Further 
details on the methodology of the color approach are in appendix I. 

4.7 SPECTRAL LINEAR DISCRIMINANT- ADJUSTED NORMAL MODEL 

The basic concept of this model is that the linear discriminant analysis of a 
set of training data can provide spectral stage estimates which can be com- 
bined with the normal model to provide improved estimates of growth stages. 
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The independent performance of the linear discriminant coefficients in pre- 
dicting coincident ground truth biostages was also measured. Further details 
on the development of these models are given in appendix J. 

4.8 COMMENTS ON ERROR ESTIMATES 

All of the adjusted-normal and spectromet models were programmed to estimate 

the Julian date on which the observed ground truth biostages would occur. The 

estimates are for the midpoints of all the stages except planting and harvest, 

when beginning and end days of the crop calendar sequence are used. Each 

ground truth development stage has a finite length ranging from less than 2 

days to nearly 13 days. If the ground truth observation is taken randomly 

from a uniform distribution, the variance, due to assuming it occurs at the 
2 

midpoint, is I /12, where I is the duration of the stage. This "rounding 
error" ranges from near 0.5 days [root mean square (RMS)] to 3.75 days, 
depending on the stage. For the ground truth evaluation data set used in this 
study, the error for an algorithm that correctly estimated the stage of each 
observation would be approximately 2 days (RMS). 

The color sequence logic and the linear discriminant models were programmed to 
estimate the Feekes growth stage in each field for each Landsat acquisition. 
The error of estimate for these models was calculated by converting both the 
estimate and the ground truth data to the percentage of normal development 
scale (see appendix 0). 
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5. RESULTS 


The normal and agromet planting date estimation models were evaluated using 
the Intensive Test Site (ITS) data set. 

The predicted planting dates using the three estimators are presented in 
table 5-1. Error statistics were calculated for each field. Table 5-2 shows 
a summary of these statistics at the segment level for all segments. 

No tests of significance have been performed. However, the Feyerherm model is 
clearly superior to the Stuff-Phi nney model and shows a sizable reduction in 
bias and root mean square error (RMSE) when compared to the normal model. 

Figure 5-1 indicates that the estimates of the Stuff-Phi nney model are 
strongly correlated with those of the normal model, whereas the Feyerherm 
estimates are relatively independent. The Stuff-Phinney model was eliminated 
from further consideration. The Feyerherm model was selected as the best 
agromet planting date model for use in the evaluation of the model systems. 

The results of the comparative evaluation of the growth stage estimation 
models are summarized in table 5-3. The table shows results for the normal, 
spectromet, agromet, and adjusted- normal models when they were applied to the 
TY data set. The results of the Badhwar-normal model are given for the fields 
upon which it could be applied. The results of direct growth stage estimates 
using the color sequence logic and the linear discriminant model are also 
shown. These models were evaluated on those observations which had coincident 
Landsat and ground-truth observations. Table 5-4 shows the performance of the 
seven models by segment and state. The remarks in the following paragraphs 
are based on these tables. 

Evaluation of planting date estimation approaches may be made through compar- 
ison of their results when used in conjunction with the normal model as a 
biostage estimator. The discriminant spectral starter algorithm performed the 
best, showing over a 50 percent reduction in variance as compared to that of 
the normal model when averaged over all segments. 
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TABLE 5-1.- PREDiaED PLANTING DATES USED IN EVALUATION 
OF NORMAL AND METEOROLOGICAL STARTER MODELS 


Segment 

Year 

Normal 

Feyerherm 

Stuff-Phinney 

1965 

1975 

143 

144 

157 

1987 

1975 

137 

136 

149 

1687 

1976 

111 

IDO 

126 

1965 

1976 

143 

126 

156 

1966 

1976 

143 

129 

155 

1967 

1976 

143 

124 

154 

1970 

1976 

125 

< > 109’ 

141 

1987 

1976 

137 

^ > 12ll' 

150 

1687 

1977' 

111 

107 

124 

1987 

1977 

137 

121 

143 


1 

f 
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TABLE 5-2.- RESULTS OF STARTER MODEL EVALUATION 


Segment 

Year 

D 

Normal 

Feyerherm 

Stuff-Phi nney 

RMSE 

Bias 

(b) 

RMSE 

Bias 

(b) 

RMSE 

Bi as 
(b) 


wss 

45 

5.993 

2.489 

5.651 

1.489 

15.339 

14.178 



8 

8.254 

.125 

8.329 

1.125 

16.359 

14.125 


1976 

5 

5.495 

/ 5 ;000 

14.184 

14.000 

16.162 

16.000 

1970 

1976 

2 

15.652 

14.000 

6.964 

-2.50 

30.806 

30.000 

1965 

1976 

8 

4.213' 

-1.250 

18.688 

-18.250 

12.420 

11.750 

1987 

1976 

12 


33.6 

18.3 

17.6 

46.9 

46.6 


1976 

31 


19.710 

10.477 

5.710 

32.904 

31.710 


1976 

9 

23.840 

19.889 

15.868 

8.889 

37.283 

34.889 

1687 

1977 


9.654 

-5.00 

12.215 

-9.00 

11.500 

8.00 

1987 

1977 


31.700 

31.182 

16.217 

15.182 

37.616 

37.182 

Total 



18.518 

11 .647 

11.900 

3.831 

27 .486 

23.619 


= number of fields used for evaluation. 

^Error defined as estimated planting date minus observed planting date. 
Positive bias indicates model is late, while negative bias indicates m^el 
is early. — ' 
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Figure 5-1.- Relationships of normal and meteorological model estimates of planting dates. 


TABLE 5-3.- SUMMARY OF STATISTICS ON RESULTS OF APPLYING BIOSTAGE 
ESTIMATION MODELS TO THE GROUND TRUTH EVALUATION DATA SET 
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(N = 464) 



TABLE 5-4.- STATISTICAL COMPARISON OF BIOSTAGE ESTIMATES WITH 
FIELD GROUND TRUTH. BY SEGMENT. STATE, AND OVERALL^ 
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Tata) nunber of observattons per state. 






The other spectral starter models showed intermediate results between those of 
the discriminant spectral starter model and the Feyerherm agromet model. It 
seems likely that the spectral techniques pick up at least part of the real 
fi el d-to- field variaton. 

The two agromet crop development models can be compared through their perform- 
ance when combined with the same starter model. The results of the Feyerherm- 
Robertson and the Feyerherm-Doraiswamy models are shown in tables 5-3 and 5-4. 
The Feyerherm agromet starter model was used as a baseline. The Doraiswamy 
model showed a slightly better performance than did the Robertson model. 

The performance of the basic adjusted-normal and agromet systems along with 
that of the normal model is given by growth stage in table 5-5. Examination 
of the normal models indicates that, although the crop was slow in development 
before the heading stage the season as a whole was near average. During head- 
ing and the early part of grain filling, the crop developed faster than normal, 
slowed during ripening, and was harvested somewhat before the normal date. 

The adjusted- normal models had the same development pattern. The agromet 
stage model picked up part of the weather- induced development variation, with 
the Doraiswamy model performing substantially better than the Rob erts on model 
in stages from heading to ripe. Appendix K contains additional details on 
evaluations of the two agromet stage estimation models. 

Contrary to what, might be expected, the (average) number of acquisitions in a 
segment had little effect on the error. This is demonstrated in table 5-6, 
which shows a comparison of the acquisition number with the RMSE, using the 
two models with the lowest RMSE overall. The explanation may be that the 
spacing of the acquisitions rather than their absolute number affects 
performance. 

The discriminant spectral starter model was combined with the Doraiswamy 
model. This spectromet model showed the best performance of those models 
evaluated in this study (tables 5-3 and 5-4). 
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TABLE 5-5.- SUMMARY OF STATISTICS ON RESULTS OF APPLYING BIOSTAGE ESTIMATION MODELS 
TO GROUND-TRUTH EVALUATION DATA SET BY GROWTH STAGE 
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^Due to computlonal problems, these statistics were calculated using the median ground truth for each segment. 




TABLE 5-6.- EFFECT OF NUMBER OF ACQUISITIONS ON THE 
RMSE OF THE DISCRIMINANT-NORMAL AND 
DISCRIMINANT-DORAISWAMY MODELS 


Number of 
acquisitions 

Number of 
segments 

RMSE, 

discriminant 
and normal 

RMSE, 

discriminant and 
Dorai swamy 

3 

1 

16.0 

9.9 

4 

5 

9.1 

8.3 

5 

13 

10.3 

10.2 

6 

8 

9.6 

9.8 

7 

1 

8.2 

8.7 
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The two spectral stage estimation models showed similar performances. A sum- 
mary of the overall results is given in table 5-3. Table 5-7 provides addi- 
tional details on the results obtained with these models at the segment and 
state levels. The results of the Badhwar multitemporal stage estimator, which 
uses the date of peak greenness to adjust the normal crop calendar, are also 
shown. Comparisons of the results of these models and the results of the 
other systems may be made only indirectly due to the different sample sizes 
available for evaluation. 

The following general observations on model performance were made. 

a. The poor performance of the Robertson model is due largely to failure to 
detect more rapid than normal crop development in most segments during the 
heading and flowering stages. 

b. The Doraiswamy model picks up some but not all of this speed-up in devel- 
opment. Future incorporation of a stress variable may improve the 
accuracy of the model. 

c. The relatively poor performance of the adjusted-normal models is due 
largely to the fact that they are coupled with the rigid normal model. 

d. However, even when the spectral models are used independently (on those 
observations which have both ground truth and spectral data) the accuracy 
is not high (approximately 11 days RMS), apparently because spectral 
growth stages correspond only roughly to physiological growth stages. 

This is to be expected since, for example, stress can produce premature 
chlorophyll disappearance while excess nitrogen can delay it. Similarly, 
different varieties could well have different spectral properties at a 
given physiological stage, or similar properties at different stages. 

e. The real strength of the spectral models has not been measured in the 
evaluation technique employed in this study. This strength lies in their 
ability to identify multitemporal spectral sequences, which analysts tend 
to refer to as spectral crop calendars. All three spectral approaches can 
capture these multitemporal patterns. 
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TABLE 5-7.- SUMMARY OF STATISTICS FOR FITS OF 
BIOSTAGE ESTIMATION MODELS 


State 

Segment 

Badhwar 

Color 

Oiscrimi nant 

N 

(4) 

Bias 

RMSE 

a 

Bias 

RMSE 

N 

(a) 

Bias 

B 

N. Oak. 

1387 

0 



5 

12.131 

-5.78 

5 


7.204 

N. Oak. 

1392 

25 

-0.80 

9.592 


10.798 

-0.80 

32 


14.466 

N. Oak. 

1394 


2.00 

7.228 

12 

5.693 

1.69 

16 


4.684 

N. Oak. 

1457 

31 

3.90 

13.109 

17 

9.720 

-0.61 

21 

1.75 

4.720 


1461 

61 

0.05 

10.063 

60 

9.900 

-2.06 

38 

-10.34 

18.850 


1472 

0 



2 

5.890 

4.16 

3 

3.48 

7.685 


1473 

16 

1.88 

4.757 

0 



4 

-3.49 

4.737 


1485 

m 



2 

2.810 

-2.81 

2 

8.99 

8.986 

iWB 

1518 

73 

7.38 

14.008 

55 

12.591 

1.90 

66 

-0.79 

9.361 

Mont. 

1537 

11 

-9.27 

11.778 

9 

14.166 

-0.35 

7 

2.14 

6.943 

Mont. 

1542 

29 

-5.79 

11.441 

17 

21.536 

-10.29 


-4.32 

19.598 

Mont. 

1544 

18 

-5.06 

10,195 

10 

7.741 

0.93 


4.86 

10.930 

Mont. 

1553 

18 

14.06 

16.487 

11 

16.686 

3.90 

9 

4.13 

6.789 

Minn. 

1566 

42 

-5.50 

12.552 

30 

7.303 

3.12 

30 

6.98 

11.922 

N. Dak. 

1584 

22 

7.59 

10.359 




4 

-1.98 

2.287 

N. Oak. 

1612 

27 

4.59 

9.907 

12 

7.400 

-2.64 

12 

-3.09 

4.846 

N. Oak. 

1619 

35 

7.31 

10.709 

16 

9.756 

-2.08 

16 

-1.65 

5.628 

N. Oak. 

1636 

49 

-4.35 

8.872 

49 

12.311 

-1.84 

40 

-1.77 

8.747 

N. Dak. 

1658 

0 



11 

9.070 

-1.63 

12 

-0.49 

1.684 

N. Oak. 

1664 

6 

4.17 

6.940 

11 

9.997 

-3.33 

12 

-0.79 

4.952 

S . Dak . 

1668 

33 

-1.76 

8.385 

18 

8.737 

1.20 

24 

8.44 

12.287 

S . Oak ^ 

1755 

0 



5 

12.926 

-1.84^ 

- 8 

0.73. 

4.883- 

S. Dak. 

1811 

54 

5.44 

8.840 

10 

9.990 

8.83 

12 

8.41 

10.613 

Minn. 

1825 

52 

-0.94 

10.613 

23 

4.770 

2.85 

26 

-2.45 

5.836 

N. Dak. 

1918 

5 

-0.20 

3.376 

4 

10.256 

-0.14 

4 

7.58 

9.426 

N. Dak. 

1920 


9.00 

14.459 

8 

16 .498 

11.51 

9 

6.55 

7.686 

N. Oak. 

1924 


0.71 

9.653 

12 

16.735 

-15.31 

12 

-2.51 

4.301 

Mont. 

1942 

■■Km 



12 

13.781 

-1.31 

12 

9.09 

14.990 

Mont. 


76 

-1.42 

12.610 

59 

2.44. 

10.020 

56 

2.52 

14.20 

N. Oak. 


351 

2.58 

10.458 

251 

-2.40 

16.256 

240 

-2.37 

10.58 

S. Oak. 


87 

2.71 

8.670 

35 

-1.84 

10.929 

46 

7.12 

10.76 

Minn. 


167 

1.55 

12.668 

108 

2.72 

9.596 

122 

89.61 

9.47 

Overall 


681 

1.90 

11.094 

453 

-0.55 

11.474 

464 

■D 

10.8 


» number of observations. 
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In other words, all the models show that most of the fields follow a similar 
spectral sequence. What is not known is the extent to which these models can 
indicate contrasting patterns for nonwheat or nonsmall grain fields. The best 
way to test the models for capability in separating land-use categories is to 
use a ground truth data set involving several categories other than wheat. 
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6. CONaUSIONS AND RECOMMENDATIONS 


The results of this evaluation are not sufficiently conclusive to warrant the 
elimination of any of the models except the Stuff-Phinney starter model. The 
lack of definitive differences is due to several factors. First, no planting 
dates were available in the data set. Second, a majority of the segments did 
not depart substantially from normal, as evidenced by the relatively small 
figures for the normal model as shown in table K-1 (see appendix K) . Third, 
the models that performed the best were those that had been fit, at least in 
part, on the same data set. 

Because of the inconclusiveness of the results, it is recommended that all 
models be carried forward for independent tests using the 1979 blind site data 
set . 

A further recommendation is that the spectral approach to planting date, esti- 
mation used by the Environmental Research Instute of Michigan (ERIM) be 
evaluated in the same way as was the Badhwar model. In the ERIM spectral 
approach, peak greenness is used as a factor. Further testing of the Badhwar 
model should be delayed until the ERIM evaluation has been completed. 

Finally, the 1979 data set should also be used to test the barley version of 
the Robertson model. - 
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Figure A-2.- Segment location map for South Dakota. 
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Figure A-3,- Segment location map for Montana. 



OmG't?5AL PAGE ^ 
OF FOC?l QUALITY 


mihm»:scta 



Figure A-4.- Segment location map for Minnesota. 
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APPENDIX B 


THE FEYERHERM STARTER MODEL 

Development of the Feyerherm spring wheat starter model (ref. 3) was based on 
the concept that a cool and wet early spring could delay planting while a warm 
and dry condition permits early planting. The weather variations among seasons 
and regions was defined as a wanning and planting (W/P) day. Each W/P day is 
assigned a 0 to 1 value for each calendar day beginning January 19 (which is 
arbitrarily chosen and coincides roughly with the coldest time of the year in 
the northern hemisphere). Accumulated W/P days were then related to the per- 
centage of wheat planted. When the accumulated W/P value reached 35.5, it was 
assumed that '50 percent of the crop had been planted. 

The general form of the model was as follows: 

W/P = 0 if TA < 32 

W/P = a (TA-32) (PRE) if 32 < TA < 32 + 1/a 
W/P = 1 if TA > 32 + 1/a 


where 

TA (°F) = the daily average air temperature [=»(max. + min.)/2] 
a - the threshold value 

PRE = a value between 0 and 1 as a function of the previous 3 days of 
precipitation 

For spring wheat, Feyerherm assigned values of a = 0.1 and PRE = 1 since 
precipitation had an effect that was statistically insignificant. 
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APPENDIX C 


THE STUFF-PHINNEY STARTER MODEL 

For the Stuff-Phinney starter model (ref. 4), the 50 percent planting date for 
spring wheat was calculated from a regression equation using daily average air 
temperature (TA), precipitation (P), and normal planting date. The daily rate 
of planting (R) was defined as follows: 

R = -0.77 + 0.045(TA) - 0.032(P) + 0.053(N) 

where N equals the actual date minus the normal planting date. The summation 
of daily rates was made only for positive values of R. 
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APPENDIX D 


STANDARDIZED CROP DEVELOPMENT SCALE 

Growth and development of spring wheat has been observed and described accord- 
ing to one of a variety of growth scales, including the Feekes scale, the 
Robertson scale, the ITS scale, and various stage descriptions used in state 
observation programs. Most of the stages in these scales can be assigned to 
one of the 24 Feekes scale stages (see table D-1 and figure D-1). 

Crop development rate may also be estimated by phenology models from meteoro- 
logical, spectral, or historical data. Because the Feekes stages vary in dur- 
ation from about 2 days (stage designation 10.5.2-10.5.3) to about 10 days 
(stage designation 10.0-10.1), a difference of 5 days in an observed stage 
date and an estimated stage date may mean that there is no error in stage 
estimation or that there is an error of several stages. Comparison of 
different phenology models and testing against observed stage data can be 
facilitated by assigning numerical values to each stage on an approximately 
linear scale. 

In figure D-1, a value for each Feekes scale stage has been derived from the 
average portion of the planting-to-harvest period required for spring wheat to 
reach that stage (i.e., planting-to-heading = 54.63 percent of the development 
period from planting to harvest). These values are based primarily on data 
from the U.S. Great Plains spring wheat region, where dates for planting, 
heading, and harvest stages (and sometimes the ripe stage) are recorded 
annually on a CRD level. 

Percentage values were calculated for each of these stages (planting = 0 per- 
cent, heading = 54.63 percent, ripe = 99.44 percent, and harvest = 100 per- 
cent) based on 5 years of CRD-level data from Montana and North Dakota and on 
about 20 years of CRD data from Minnesota and South Dakota. Other stages were 
placed on the scale according to scattered historical CRD observations, exper- 
iment plot data, and by interpolation. The Robertson and ITS stages were 
assigned Feekes scale analogues based on definitions of the stages. 
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TABLE D-1.- RELATIONSHIPS BETIaIEEN FEEKES STAGES, ITS STAGES, AND 
ROBERTSON STAGES, WITH ESTIMATES OF PERCENTAGE OF NORMAL 
GROWTH SEASON OF EACH STAGE FOR SPRING WHEAT 


Normal percentage of season 

Feekes 

stages 

ITS 

stages 

Robertson 

stages 


Stage 

beginning. 

Stage 

midpoints 

Description 

0.0 

0.0 

0 

2 

1.0 

PI anted 

8.33 

13.88 

1 

3 

2.0 

Emergence 

19.44 

20.83 

2 

- 

2.2 

Beginning tillering 

22.22 

23.61 

3 

4 

2.4 

Till ers formed 

25.0 

26.85 

4 

- 

2.6 

Beginning pseudostem 
erection 

28.7 

30.55 

5 

- 

2.8 

Pseudostem strongly 
erected 

32.41 

33.80 

6 

- 

3.0 

Jointing 

35.19 

37.04 

7 

• 

3.1 

Second node formed; 
second-to-last leaf 
visible 

38.89 

40.28 

8 

- 

3.25 

Last leaf visible 

41.67 

43.52 

9 

- 

3.35 

Last leaf ligule 
visible 

45.37 

50.00 

10.0 

5 

3.5 

Boot 

54.63 

56,02 

10.1 

6 

4.0 

First heads just 
visible 

57.41 

58 .-80 


■- 

4.1 

One-fourth headed 

60.19 

61.11 


- 

4.2 

One- half headed 

62.04 

63.42 

10.4 

7 

4.25 

Three-fourths headed 

64.81 

71.30 

10.5 

- 

4,35 

All plants headed 

77.8 

80.11 

11.1 

- 

4.8 

Milky ripe 

82.41 

84.70 

11.2 

8 

5.0 

Soft dough mealy ripe; 
kernel soft but dry 

87.0 

90.72 

11.3 


5.4 

Kernel hard, but 
difficult to divide 
with thumbnail 

94 .44 

97.22 

11.4 

9 

6.0 

Ripe 

100.0 

100.0 

11.5 

10 

7.0 

Harvested 
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Figure D-1.- Spring wheat growth scales conversion chart. 
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TABLE E-1.- NORMAL SPRING WHEAT PLANTING AND HARVEST 
DATES FOR SEGMENTS IN EVALUATION DATA SET 


Segment 

State 

Planting 

Harvest 

Duration 

1387 

N. Dak. 

140 

235 

95 

1392 

N. Dak. 

141 

240 

99 

1394 

N. Dak. 

143 

242 

99 

1457 

N. Dak. 

143 

242 

99 

1461 

N. Dak. 

141 

240 

99 

1472 

N. Dak. 

134 

228 

94 

1473 

N. Dak. 

134 

228 

94 

1485 

S. Dak. 

117 

219 

102 

1518 

Minn. 

137 

243 

106 

1537 

Mont. 

128 

236 

108 

1542 

Mont , 

128 

236 

108 

1544 

Mont. 

128 

236 

108 

1553 

Mont. 

121 

231 

no 

1566 

Mi nn . 

119 

224 

105 

1584 

N. Dak. 

140 

235 

95 

1612 

N. Dak. 

141 

240 

99 

1619 

N. Dak. 

140 

235 

95 

1636 

N. Dak. 

135 

227 

92 

1658 

N. Dak. 

122 

219 

97 

1664 

N. Dak. 

122 

219 

97 

1668 

S. Dak. 

117 

219 

102 

1755 

S. Dak. 

111 

207 

96 

1811 

S. Dak. 

no 

209 

99 

1825 

Mi nn . 

137 

243 

106 

1918 

N. Dak. 

127 

222 

95 

1920 

N. Dak. 

127 

222 

95 

1924 

N. Dak. 

122 

219 

97 

1942 

Mont. 

128 

236 

108 
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APPENDIX F 


ROBERTSON MODEL 

Robertson's biometeorological time scale model (ref. 5) is based on 
photothermal effects on crop phonological processes. The nonlinear responses 
of photoperiod and day and night temperatures are represented in a tri qua- 
dratic format. The effects of each environmental parameter are considered 
over short phonological periods (stages) when the physiological processes are 
uniform. These stages are planting, emergence, jointing, heading, soft dough, 
and ripe. Intermediate stages corresponding to those presented in the ground 
truth (Feekes scale) were also computed by interpolation from the normal crop 
calendar. 

Maturity from one stage (SI) to the next (S2) is computed in the model as 
follows; 

12 

1 = Maturity = 22 
SI 

' '’o’ * j 

where P is the daylength, T^^gj^ is the daily maximum temperature, and is 
the daily minimum temperature. 

Coefficients ag, a^ a 2 » bg, b^, b 2 , and d 2 are determined by' an iterative 
regression technique which provides the best relationship between the three 
environmental factors and their interactions for the data set used. The 
coeffients are tabulated in table F-1. 

The model was developed for Canadian spring wheat, and there are some 
weaknesses in the model that may limit its predictive capability in other 
locations. The responses of the model to the three environmental variables 
are not similar to expected responses for spring wheat varieties grown in the 
Midwest. One evidence of this is the prediction of almost the same duration 
from planting to ripe for a range of planting dates. 
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TABLE F-1.- REGRESSION COEFFICIENTS FOR THE ROBERTSON MODEL 
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DORAISWAMY MODEL 

Examination of intermediate development stages predicted by the Robertson 
model and the unrealistic responses produced by the input parameters during a 
given crop stage indicated a need for modification of the model. The 
Robertson model computes the development period by combining all of the 
parameter responses, irrespective of their physiological or phenol ogical 
significance. For example, the photoperiod effect on phenology is carried 
throughout all stages of plant development, although it has an insignificant 
effect after flowering. In general, the predictions of the Robertson model 
for the entire development period are close to the ground truth observations; 
however, the intermediate stage predictions are sometimes quite erroneous. 

The Doraiswamy model, while using the same computational form of the Robertson 
model, separates parameter effects during each development stage. Certain 
coefficients are modified by comparison with the ground truth for a few ran- 
domly selected segments. Consequently, the new coefficients give more 
relevant weighting to the parameters; e.g., photoperiod is neglected (a^ = 0 
and Si 2 = 0) during periods in which the crop is insensitive to it. The coef- 
ficients are given in table G-1. 



TABLE G-1.- COEFFICIENTS FOR THE DORAISWAMY MODEL 
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APPENDIX H 


BADHWAR-ADJUSTED NORMAL MODEL 

The technique involved in the Badhwar-adjusted normal model consists of 
fitting field averages of Kauth greenness over at least three acquisition 
dates using the nonlinear expression 

p(t) = A(t/tQ)“ exp {0tQ [1 - (t/tg)^]} 

where p is the spectral variable (Kauth greenness), t is time (Julian date of 
acquisition). A, a, e, and tg are coefficients to be fit, and exp denotes the 
exponential function. The acquisitions must be in the postemergence to pre- 
harvest period for an acceptable fit to the assumed function. Initial values 
for the coefficients and an acceptable step size which provide good conver- 
gence for the iterative fitting procedure are also needed. These have been 
chosen to give good results over a wide range of segments, and the values used 
for the evaluation data set are given in table H-1. 

TABLE H-1.- BADHWAR INITIAL CONDITIONS 

Values Statistics 

Initial coefficients a = 18.0, B = 2.0, tg = 1.5 

■ Coefficient step change 1. x lO"® - - - - 

2 ti 

X convergence difference 1. x 10"^ 
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APPENDIX I 


COLOR-ADJUSTED NORMAL MODEL 


The first step in the color adjusted normal model (ref. 7) normalizes Landsat 
pixel channel values to the segment channel means. Initially, one divides all 
pixel (or field mean) channel values (S^) by scene (segment) channel means 
(CH^O and multiplies them by five to obtain relative energies (X^-)« 


where i = 1 for channel 1, i =2 for channel 2, and i = 3 for channel 4. 

The resulting relative energies are transformed to the three color components 
of hue, value, and chroma as follows: 


A =^X + 2 X + 3 X 
S = X^ + X 2 + X 3 


Hue = COS 


-1 


2X3 - X^ 


2VA - X^X 2 - X^X 3 - X 2 X 3 


Value = S/3 
Chroma = Va - 1/3S^ 


Pixels (or fields) are then classified into the generalized color groups 
defined in table I-l. Additional groups may be defined as needed. For each 
field and spectral acquisition, the percentage of each color is calculated and 
subjected to a decision logic based on the temporal distributions of the per- 
centages to estimate either the corresponding physiological stage midpoint or 
the biowindow equivalent. 


I-l 



TABLE I-l.- COLOR COOES DETERMINED FROM HUE, 
VALUE, AND CHROMA 


Code 

Hue in degrees 

Value 

Chroma 

Green (G) 

-120 

to 

-180 

0 to 

10 

0-8.165 

Blue (B) 

120 

to 

180 

0 to 

10 

0-8.165 

Purple (P) 

60 

to 

120 

0 to 

10 

0-8.165 

Light Red (L) 

0 

to 

60 

5 to 

10 

0-8.165 

Dark Red (R) 

0 

to 

60 

LT 5 


0-8.165 

Orange (0) 

0 

to 

-60 

0 to 

10 

0-8.165 

Yellow (Y) 

-60 

to 

-120 

0 to 

10 

0-8.165 


The stage estimate is extrapolated backward with the normal model to a 
planting date. The extrapolated values are averaged to predict a planting 
date which adjusts the normal model to predict occurrences of ground truth 
stages. 

The biowindow estimation logic is given in table 1-2, while table 1-3 shows 
the stage estimation logic. The color codes and stage codes used in both 
these tables are fully identified and described below table 1-2. 
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TABLE 1-2.- LOGIC FOR SEPARATION OF SPRING SMALL GRAIN SPECTRAL 
RESPONSES INTO BIOWINDOW EQUIVALENTS USING THE COLOR MODEL^ 
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TABLE 1-3.- LOGIC FOR SEPARATION OF SPRING SMALL GRAIN SPECTRAL RESPONSES 

INTO STAGES USING THE COLOR MODEL® 
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percentage of all pixels that are green Late - 75.0% on normal scale 

percentage of all pixels that are blue . DR-RH - 87.5% on normal scale 

percentage of all pixels that are orange RH - 100.0% on normal scale 

percentage of all pixels that are yellow 
percentage of all pixels that are orange or yellow 
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APPENDIX J 


LINEAR DISCRIMINANT BIOSTAGE ESTIMATION 

The potential for use of linear discriminant functions in identification of 
crop development stages was demonstrated for wheat by Phinney et al. (ref. 8). 
For application to spring wheat, a training set consisting of coincident 
ground data and spectral data was compiled. The data were collected at 35 
spring wheat segments in Montana, North Dakota, South Dakota, and Minnesota 
during the 1978 crop year. The total number of observations was 590. 

Table J-1 shows the number of observations available for each observed 
development stage. Data were pooled for Feekes stages 10.2 to 10.4. Although 
the number of available samples was small for several stages, the effect of 
pooling additional stages was not investigated. Eighteen different stage 
classes were considered. 

Six different spectral transforms were used in addition to Sun-angle correction 
of the raw Landsat multispectral scanner (MSS) channel data. These transforms 
included the tassel led cap greenness and brightness, the normalized difference, 
the perpendicular vegetation index (PVI), relative energy, and the cubic 
colors. The spectral variables were summarized at the field level for each 
acquisition date. 

Preliminary studies showed that fields near the beginning and end of the 
growth cycle were often confused when only spectral data were used. Hence, 
the difference between the acquisition date and the normal planting date for 
the segment's CRD provided a segment-level variable to complement the field- 
level variable. 

The results of the fit to the training data set are shown in table J-2. The 
percent correct classification for each set of variables is shown for Feekes 
stage and by biowindow. Table J-3 shows the magnitude of the misclassifica- 
tion error in days. The normal growth season was used to convert error in 
stages to error in days. 


J-1 



TABLE J-1.- SUMMARY OF OBSERVATIONS USED FOR FITTING OF 
LINEAR DISCRIMINANT COEFFICIENTS 
BY DEVELOPMENT STAGE 


Feekes stage 


Number of observations 



^These three observations were pooled for Feekes stages 
10.2, 10.3, and 10.4. 
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TABLE J-2.- RESULTS OF LINEAR DISCRIMINANT FITS TO THE TRAINING DATA 


[Based on 590 observations] 


Spectral 

variable 

Number of 
variables 

Percent correct 

classification 

Feekes 

stage 

Bi owi ndow^ 

1 

2 

3 

D 

All 

Color 

8 

53.9 

93.1 


63.4 

93.6 

82.9 

MSS channel 

5 

54.7 

91.9 


73.9 

93.6 

84.6 

Relative energy 

5 

54.1 

94.2 

74.8 

73.2 

90.0 

83.7 

Tasselled cap 

3 

53.6 

90.8 

74.1 

67.6 

94.3 

82.2 

Normalized difference 

2 

51.2 

91.9 

74.8 

60.6 

93.6 

75.8 

PVI 

2 

50.8 

84.4 

77.8 

61.3 

91.4 

79.0 


^Biowindow definition: 

1 * Feekes stage 1-2 

2 s Feekes stage 3-10.5 

3 s Feekes stage 11.1-11.3 

4 « Feekes stage 11.4-11.5 

TABLE J-3.- MISCLASSIFICATION STATISTICS IN 
DAYS FOR RESULTS OF LINEAR DISCRIMINANT 
FITS TO THE TRAINING DATA 

[Based on 590 observations] 


Spectral variable 

Bias 

RMSE 

Color 

-1.1 

12.41 

MSS channel 

-.9 

9.72 

Relative energy 

-.9 

9.46 

Tasselled cap 

-.6 

9.79 

Normalized difference 

-1.0 

10.90 

PVI 

-.2 

12.97 
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The coefficients obtained from the training data were applied to a larger,^ 
ground truth evaluation data set with 1057 observations. The training- ‘data 
set was a subset of the larger evaluation data set. There were 464 coincident 
observations of spectral ground data. Table J-4 shows the performance of the 
candidate model on the coincident data. 

Estimates of planting dates were made. for each field in each spectral acquisi- 
tion. The estimated stage was calculated from the estimated stage by adjust- 
ing the normal crop calendar. A final planting date estimate for each field 
was made by averaging the individual acquisition estimates. 

It should be noted that the date that a given ground truth stage occurred can 
be estimated by adjusting the normal crop calendar by the estimated planting 
date. The results of this approach are also shown in table J-4. Comparison 
of the direct estimates and the adjusted normal estimates indicate that no 
skill is lost by using the spectral model solely to estimate planting date. 

The linear discriminant models were used to estimate planting dates for each 
field using all available spectral acquisitions in the ground truth evaluation 
data set. The estimated planting dates adjusted the normal model to estimate 
the date of occurrence of ground truth stages in the data set. There were 
1057 such ground observations. Table J-5 shows a summary of the results. 

Reviewing tables J-3, J-4, and J-5, it may be seen that the relative energy 
spectral transform provided the best results, and the tassel led cap and cubic 
color transform performed nearly as well. The relative energy transform was 
selected for comparison with other techniques, and its linear discriminant 
coefficients are given in table J-6. 

For the relative energy transform, on each acquisition date, the variables for 
each field are calculated using the following methodology: 

X(l) = relative energy for MSS channel 1 

X(2) = relative energy for MSS channel 2 

X(3) = relative energy for MSS channel 3 

X(4) = relative energy for MSS channel 4 

X(5) = acquisition date - normal planting date 
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TABLE J-4.- MISCLASSIFICATION STATISTICS IN DAYS FOR RESULTS 
OF LINEAR DISCRIMINANT MODELS APPLIED TO 
GROUND TRUTH EVALUATION DATA 


[For 464 observations] 


Spectral variable 

Di rect 

Estimated planting 
date + normal 

Bias 

RMSE 

Bias 

RMSE 

Color 



0.3 

11.05 

MSS channel 

B 


.2 

10.43 

Relative energy 

.0 


.1 

10.08 

Tassel led cap 

B 

10.83 

-.2 

10.27 

Normalized difference 

M 

11.86 

1.0 

11.39 

PVI 

H 

14.05 

.3 

12.73 


TABLE J-5.- SUMMARY OF STATISTICS FOR ADJUSTING NORMAL 
MODELS WITH PLANTING DATE ESTIMATES DERIVED 
FROM THE LINEAR DISCRIMINANT MODELS 


Spectral variable 

Bi as 

RMSE 

Color 

-0.9 

10.85 

MSS channel 

-0.1 

11.97 

Relative energy 

-1.5 

10^32 

Tasselled cap 

-1.1 

10.67 

Normalized difference 

0.2 

12.54 

PVI 

0.1 

13.18 
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TABLE J-6.- LINEAR DISCRIMINANT COEFFICIENTS FOR RELATIVE ENERGY 


Group 

Feekes 

stage 

Relative energy 

by MSS channel 

Days since 
planting 


1 

2 

3 

4 

Constant 

1 

1.0 

80.42 

-30.00 

-54.79 

43.35 

0.08 


2 

2.0 

83.65 

-32.43 

-52.22 

41.56 

0.16 

-107.80 

3 

3.0 

82.60 

-31.94 

-53.50 

43.42 

0.19 

-108.44 

4 

4.0 

82.77 

-33.41 

-49.65 

39.97 

0.28 

-107.15 

5 

5.0 

82.54 

-34.89 

-46.30 

38.21 

0.25 

-106.68 

6 

6.0 

78.91 

-31.36 

-52.97 

44.28 

0.27 

-102.98 

7 

7.0 

80.29 

-33.63 

-49.97 

43.45 

0.29 

-112.92 

8 

8.0 

79.69 

-33.40 

-52.92 

48.17 

0.34 

-126.25 

9 

9,0 

77.61 

-31.28 

-58.27 

52.64 

0.38' 

-122.92 

10 

10.0 

81.12 

-34.02 

-48.64 

41.76 

0.38 

-117.10 

U 

10.1 

76.50 

-30.93 

-52.50 

45.43 

0.35 

-106.28 

12 

10.2-10.4 

79.46 

-31.85 

-53.37 

44.05 

0.34 

-103.03 

13 

10.5 

77.06 

-30.25 

-59.11 

50.74 

0.43 

-110.76 

14 

11.1 

74.61 

-28.19 

-62.16 

53.19 

0.48 

-108.96 

15 

11.2 

74.74 

-27.83 

-58.40 

48.45 

0.50 

-107.02 

16 

11.3 

70.91 

-25,86 

-54.86 

45.42 

0.51 

-102.41 

17 

11.4 

70.19 

-25.35 

-54.32 

45.09 

0.59 

-109.48 

18 

11.5 

76.84 

-27.91 

-56.11 

45.65 

0.68 

-133.34 
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Calculate 


S(i) = Const(i.) + X(j)Coeff (i ,j ) 

where Const(i) is the constant for the 1th group and Coeff(i,j) is the coef- 
ficient for the ith group and jth variable. 

The group which has the maximum likelihood of containing the observation will 
have the maximum value of S(i). Thus, the Feekes growth stage associated with 
the group with the maximum S(i) is the stage estimate. 

The planting date is estimated as follows: 

planting date = acquisition date - NP x DUR 

where NP is the percentage of the normal growing season to the midpoint of the 
estimated stage and DUR is the normal length of the growing season for the CRD 
containing the field. 

The average planting date obtained from all acquisitions is then the final 
estimate of the planting date. In general, acquisitions outside the normal 
growing season should not be used. 
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APPENDIX K 


EVALUATION OF THE ROBERTSON AND DORAISWAMY STAGE ESTIMATION MODELS 

The performances of the Robertson and Doraiswamy models are briefly discussed 
in conclusions 1 and 2 of the results section of this document. For a better 
understanding, the predictions of the two models and the normal crop calendar 
were tabulated for 25 days before and after the normal planting date. 

Tables K-1, K-2, and K-3 illustrate the observed improvement of the Doraiswamy 
model in segment 1918. The dates circled indicate the median occurrence of a 
particular growth stage observed in the ground truth. A crop phenology model 
having accurate predictions relative to the ground truth would have all the 
ground-truth observations distributed horizontally. Both the Doraiswamy and 
Robertson models and the normal crop calendar predictions suggest (1) a gen- 
eral increase in the development rate for crops during stages from heading to 
flowering and (2) a decrease in development during subsequent stages. The 
models are unable to effectively simulate this effect although the Doraiswamy 
model does better than both the Robertson model and the normal crop calendar. 
One of the possible explanations for this anomaly is that in 1978 a larger 
plant water stress affected the normal development of the crop. 
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TABLE K-1.- DATES OF STAGE OCCURRENCES FOR SEGMENT 1918 (GRANT COUNTY, N. DAK.) 
AS PREDICTED BY THE NORMAL MODEL FOR PLANTING DATES 
WITHIN 25 DAYS OF NORMAL 
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TABLE K-3.- DATES OF STAGE OCCURRENCES IN 1978 FOR SEGMENT 1918 (GRANT COUNTY, N. DAK.) 
AS PREDIQED BY THE DORAISWAMY MODEL FOR PLANTING DATES 
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Several of the projects in the Agricultural and Resources Inventory Surveys Through Aerospace 
Remote Sensing (AgRISTARS) program require crop phenology information, including classification, 
acreage and yield estimation, and detection of episodal events. This study evaluates several 
crop calendar estimation techniques for their potential use in the program. The techniques, 
although generic in approach, were developed and tested on spring wheat data collected in 1978. 
There are three basic approaches to crop stage estimation; historical averages for an area 
(normal crop calendars), agrometeorological modeling of known crop-weather relationships 
[agrometeorological (agromet) crop calendars], and interpretation of spectral signatures 
(spectral crop calendars). Normals serve as the baseline against which the skill of the other 
models may be tested. Agromet crop stage models require planting dates to initiate the accumu- 
lation of weather variables. When a data set does not include planting dates, these dates may 
be estimated with normal, agromet, or spectral starter models and used to start the agromet 
stage models. The 1978 spring wheat data set does not include planting dates. Consequently, 
five starter models and the normal planting date were used, including the Feyerherm and Stuff- 
Phinney agromet starter models and the Badhwar profile, color, and linear-discriminant spectral 
starter models. The planting date models were compared by using them to adjust the normal crop 
calendar to predict stages on dates for which ground truth is available. Two agromet models, 
the Robertson and Doraiswamy, were tested using normal, Feyerherm (the best agromet starter), 
and linear discriminant (the best spectral starter) planting dates. In all, 10 combinations of 
planting and biostage estimation models were evaluated. Dates of stage occurrence are esti- 
mated with biases between -4 and +4 days while root mean square errors range from 10 to 15 days. 
Results are inconclusive as to the superiority of any of the models and further evaluation of 
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